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We consider higher-order QCD corrections for vector boson production at hadron colliders. 
We present recent results on transverse-momentum resummation for Z production. Moreover 
we show numerical results from a new fully exclusive next-to-next-to-leading order (NNLO) 
calculation. 



1 Introduction 



The study of vector boson production in hadron cohisions, the well know Drell-Yan (DY) process, 
is nowadays extremely important. Owing to the large production rates and clean experimen- 
tal signatures of W and Z productions, these processes are standard candles for calibration 
purposes, they lead to precise determinations of vector boson masses and widths and provide 
important information on parton distribution functions (pdf). It is therefore essential to have 
accurate theoretical predictions for vector boson cross section and distributions. This requires 
the computation of higher-order QCD radiative corrections for such processes. 

In these proceedings we present two recent results on higher-order corrections for the DY 
process: a study of transverse-momentum (qt) resummation for Z production at the Tevatron^ 
and a fully exclusive calculation^ based on the NNLO subtraction formalism of Ref.^. 



2 Transverse-momentum distribution: fixed-order and resummation 

We consider the inclusive hard-scattering process 

hiipi) + h2{p2) ^ V{M,qT)+X ^ k+l^+X, (1) 

where hi and /12 are the colliding hadrons with momenta pi and p2, is a vector boson (which 
decays in the lepton pairs /i, ^2) with invariant mass M and transverse- momentum qx and X is 
an arbitrary and undetected final state. 

According to the QCD factorization theorem the qt differential cross section da^ /dq^ can 
be written as 




where fa/h{x,fip) are the parton densities of the colliding hadrons at the factorization scale fip, 
da^f^/dq^ are the perturbative QCD computable partonic cross sections, s (s = xiX2s) is the 
hadronic (partonic) centre-of-mass energy, and hr is the renormalization scale. 



In the region where qx ~ ?7iy, my being the mass of the vector boson {my = rawjUiz), the 
QCD perturbative series is controlled by a small expansion parameter, as{mv), and fixed-order 
calculations are theoretically justified. In this region, the QCD radiative corrections are known 
up to next-to-leading order (NLO)SI. From Fig. [T] we see that the NLO result^ is in agreement 
with the experimental data^over a wide region of transverse momenta {qx ~ 20 GeV). 

In the small-gr region (qx ^ mv), the LO and NLO calculations do not describe the data. 
This is not unexpected since in the small-gr region the convergence of the fixed-order perturba- 
tive expansion is spoiled by the presence of powers of large logarithmic terms, a'g In"^ (my / q'^) . 
To obtain reliable predictions these terms have to be resummed to all orders. 
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Figure 1: The gT-spectrum of the Drell-Yan e'^e~ pairs produced in pp colhsions at the Tevatron Run I. Theoretical 
results are shown at LO and NLO including the following scale variations: mz/2 < /iF,A*-R < 2 mz , with the 

constraint 1/2 < hf/i-Ir < 2. 




The resummation is performed at the level of the partonic cross section, which is decomposed 
as da^^/dq^ = da^^^^^'^'''^ / dq^ + dd^^^^'^'^ /dq^ . The term dcr^^^^'^^ '' contains all the logarithmically 
enhanced contributions (at small qx), which have to be resummed to all orders in as, while the 
term do^^^^^'^'^ is free of such contributions and can be evaluated at fixed order in perturbation 
theory. Using the Bessel transformation between the conjugate variables qx and 6, the resummed 
component da^f^^'^^^ '' can be expressed as 

°% {qT,M,s,as)=crlo{M)— / db - UhqT)Wl{h,M,s,as) , (2) 
aqrp s Jq z 

where a^o is the Born partonic cross section and Jo{x) is the 0-order Bessel function. By taking 
the A^- moments of W with respect to the variable z = M'^Js at fixed M, the resummation 
structure of jy can be organized in the exponential formQ 

W];;(6, M, as) = nl {as) x exp{^;v(a5, L)} , with L = IniQH^bl), bo = 26"^^ . (3) 

The scale Q ^ M ^ my, that appears in the above formula is named resummation scale and 
it parameterizes the arbitrariness in the resummation procedure. Variations of Q around my 
can be used to estimate the size of higher-order logarithmic contributions that are not explicitly 
resummed in a given calculation. 

The process dependent function Ti^ includes all the perturbative terms that behave as 
constants as qx 0. It can thus be expanded in powers of as = as{lJ^\)'- 

n-iV ( \ Ti , OLS njV {I) ( Ois\'^ njV {2) , 1 ,.s 
= [1 + — ^— j Hj^'-' + .-.y (4) 



"For the sake of simplicity, here we consider only the case of the diagonal terms in the flavour space. For the 
general case and a detailed discussion of the resummation formalism see Ref.^. 



The universal exponent Gn contains all the terms that order-by-order in as are logarithmically 
divergent as 6 — > oo (i.e. — > 0). The logarithmic expansion of reads 



GNias,L) = Lg 



{asL) + <7^') (asL) + (asL) + . . . 



(5) 



where the term L g^^^ collects the leading logarithmic (LL) contributions, the function g^^^ 
includes the next-to-leading leading logarithmic (NLL) contributions and so forth 0. 

Finally the finite component has to be evaluated starting from the usual fixed-order perturba- 
tive truncation of the partonic cross section and subtracting the expansion of the resummed part 

at the same perturbative order: [da^i^^^^'^ /dq'^]f,o. = [d^ab/^lT^f-o- ~ [d^ab^^'^^^ /^lT]f-o-- This 
matching procedure between resummed and finite contributions guarantees to achieve uniform 
theoretical accuracy over the entire range of transverse momenta. 

The inclusion of the functions g^^^ , g^^^ , "H^^^^ in the resummed component and of the finite 
component at LO (i.e. 0{as)) allows us to perform the resummation at NLL+LO accuracy. 

(3) V(2) 

The inclusion of the functions gj^ and and of the finite component at NLO leads to a full 

V(2) 2 

NNLL+NLO accuracy. Since the coefficient has been computed only recently^, here we 

limit ourselves to presenting results up to NLL+LO accuracy. 

In Fig. [2] we compare our NLL+LO resummed spectrum^ll (with different values of the facto- 
rization, renormalization and resummation scale) with the Tevatron data. We find that the 
scale uncertainty is about ±12 — 15% from the region of the peak up to the intermediate qt 
region (qt ~ 20 GeV), and it is dominated by the resummation-scale uncertainty. Taking into 
account the scale uncertainty, we see that the resummed curve agrees reasonably well with the 
experimental points. We expect a sensible reduction of the scale dependence once the complete 
NNLL+NLO calculation is available. 

We note that in Fig. [2] the theoretical results are obtained in a pure perturbative framework, 
without introducing any models of non-perturbative contributions. These contributions can be 
relevant in the qt region below the peak. 
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Figure 2: The gT-spectrum of the DreU-Yan e^e~ pairs produced in pp colhsions at the Tevatron Run I. 
Theoretical results are shown at NLL+LO, including scale variations. Left side: mz/2 < fiF,iJ.R < 2 mz , with 
the constraint 1/2 < hf/i-ir < 2. Right side: mz/4 < Q < mz 



''To reduce the impact of unjustified higher-order contributions in the large-gy region, the logarithmic variable 
L in Eq. (|3]), which diverges for b — + 0, is replaced by L = In (Q^b^/feg + l) ■ As a consequence of this replacement, 
integrating the qr distribution over we obtain the corresponding total cross section: dq^{da / dq^) nll+lo = 



3 Fully exclusive NNLO Drell-Yan calculation 



We now consider QCD radiative corrections at the fully exclusive level for the process in Eq. [TJ 
The purpose is to compute observables da^ , with arbitrary (though infrared safe) kinemati- 
cal cuts on the final-state leptons and the associated jet activity. Provided the observable is 
sufficiently inclusive over the small-g7^ region, it can reliably be computed at fixed order in 
perturbation theory. 

Following Ref.'^, we observe that, at LO, the transverse momentum qt of V is exactly zero. 
This means that if 7^ the (N)NLO contributions is given by the (N)LO contribution to 
the final state V + jei(s): daJj^^j^j^Q\q^^o = c^o'^^^^q ■ We compute da^'^^Q^^ by using the 
subtraction method at NLO and we treat the remaining NNLO singularities at qt = by 
the additional subtraction of a counter-term | constructed by exploiting the universality of the 
logarithmically-enhanced contributions to the qt distribution (see Eq. [3]) 
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where 'H^j^^j^^q is the process dependent coefficient function of Eq. [H 

We have encoded our NNLO computation in a parton level Monte Carlo event generator. The 
calculation includes finite-width effects, the 7 — Z interference, the leptonic decay of the vector 
bosons and the corresponding spin correlations. Our numerical code is particularly suitable 
for the computation of distributions in the form of bin histograms, as shown the illustrative 
numerical results presented in Fig. [3l 
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Figure 3: Left side: transverse mass distribution for W production at the Tevatron. Right side: distributions in 

PTmin and PTmax for the Z signal at the Tevatron. 
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